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An alternative theory of gravity has recently been proposed by Bekenstein, named Tensor- Vector- 
Scalar (TeVeS) theory, which can explain many galactic and cosmological observations without the 
need for dark matter. Whilst this theory passes basic solar system tests, and has been scrutinized 
with considerable detail in other weak-field regimes, comparatively little has been done in the strong- 
field limit of the theory. In this article, with Cowling approximation, we examine the oscillation 
spectra of neutron stars in TeVeS. As a result, we find that the frequencies of fundamental modes 
in TeVeS could become lager than those expected in general relativity, while the dependence of 
^\ \ frequency of higher overtone on gravitational theory is stronger than that of lower modes. These 

■ imprints of TeVeS make it possible to distinguish the gravitational theory in strong-field regime via 

' the observations of gravitational waves, which can provide unique confirmation of the existence of 

\ scalar field. 

, PACS numbers: 04.40.Dg, 04.50.Kd, 04.80.Cc, 97.60.Jd 



I. INTRODUCTION 

Tests of gravitational theories in the strong-field regime are extremely important because, unlike the weak-field, 
r : they are still largely unconstrained by observations. However, with developing technology, it is becoming possible to 
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observe compact objects with high accuracy. For example, observations of emitted X-rays and 7-rays from compact 
objects could be used to directly test the strong-field regime of a gravity theory 1]. Furthermore, future observational 
developments, for example with gravitational waves, will allow us to obtain different physical properties for compact 
objects, which will further allow for testing in the strong-field regime. From a theoretical point of view, there are 
attempts to test theories of gravity in the strong-field regime by using surface atomic line redshifts [3| or gravitational 
' waves from the neutron stars [3]. In these investigations, the possibility of distinguishing Scalar- Tensor (ST) theory, 
. proposed by Damour & Esposito-Farese Q, from General Relativity (GR) was discussed (see Psaltis [l| for a review). 
fSJ ' Whilst the existence of scalar fields has not been experimentally verified, several experiments in the weak-field limit 
. of GR set severe limits on the existence and strength of such fields j5j| . 

' Recently, Tensor- Vector-Scalar (TeVeS) theory has attracted considerable attention as an alternative gravitational 
^\ ' theory. TeVeS was proposed by Bekenstein as a relativistic theory for Modified Newtonian Dynamics Q- As 
, such, it explains galaxy rotation curves and the TuUy-Fisher law without the existence of dark matter. TeVeS has 
also successfully explained strong gravitational lensing as well as key features of the cosmic microwave background 
and galaxy distributions through an evolving Universe fiol w ithout cold dark matter. While in the strong-field 
^ ' regime of TeVeS, Giannios found the Schwarzschild solution and Sagi and Bekenstein generalized this to the 
^ ', Reissner-Nordstrom solution Furthermore, Contaldi, Wiseman & Withers have found vacuum solutions for a 

constant scalar field [3]. More recently, the Tolman-Oppenheimer-Volkoff (TOV) equations in TcVcS arc derived 



by Lasky, Sotani & Giannios with which one can produce static, spherically symmetric neutron stars, and they 
showed the possibility of distinguishing TeVeS from GR by way of redshift observations. In this article, we examine 
whether observations of gravitational waves associated with the neutron star oscillations can provide an alternative 
way of probing the gravitational theory in the strong-field regime. 

The attempt to estimate the stellar parameters, such as mass, radius and equation of state (EOS), via their 
oscillation properties is not a new idea. Helioseismology is established fields in astronomy and one could know the 
information about the interior of our sun. In the late '90s, it was suggested the possibility to reveal the compact star 
properties through the oscillation spectra which is called as gravitational wave asteroseismology. The stellar mass, 
radius and EOS can be deduced by an ana.lysis of the oscillation spectrum of fundamental, pressure and spacetime 



modes, i.e., /, p and w modes (e.g., [l^, [T? 
the examination of r mode oscillations (e.g. 



While the rotation period of a compact star can be revealed by 
3 H^), where such frequencies are proportional to the rotation rate. 
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Furthermore, the detailed analysis of the gravitational waves makes it possible to determine the radius of accretion 
disk around supermassive black hole |20| or to know the magnetic effect during the stellar collapse [2l| . 

In general, the oscillations of a neutron star in TeVeS could produce not only gravitational but also scalar and 
vector waves, which is similar to the case in the scalar-tensor theory ^22*1, and the direct detection of scalar and/or 
vector waves would be a unique probe for the gravitational theory. Still, we will show that it is not necessary to detect 
these waves, because the obvious imprints due to the existence of scalar and vector fields will be apparent in the 
spectrum of gravitational waves associated with the stellar oscillations. Although we adopt Cowling approximation 
in this article, the more complicated analysis including the metric, vector and scalar perturbations will be seen near 
future somewhere. 

This article is organized as follows. In the next section, we describe the fundamental parts of TeVeS and TOV 
equations in TeVeS, where we also show the neutron star models. In Sec. IIIII we derive the perturbation equations 
with Cowling approximation. Then the oscillation spectra of neutron stars in TeVeS are shown in Sec. IIVI finally 
we discuss the results related to gravitational wave asteroseismology in Sec. |Vl In this article, we adopt the unit 
of c G = 1, where c and G denote the speed of light and the gravitational constant, respectively, and the metric 
signature is (-, +, +, +). 



II. STELLAR MODELS IN TEVES 



TeVeS 



Since details of TeVeS can be found in 0, we only mention here the fundamental parts of the theory that are 
necessary for the present calculations. TeVeS is based on three dynamical gravitational fields; an Einstein metric g^jy, 
a timelikc 4- vector field Z^^, and a scalar field There is also a nondynamical scalar field, a. The vector field fulfills 
the normalization condition, g^JA^W^ = — 1, and the physical metric is given by 



9^ll> = e '^'-^g^.u - ^U^Ui, sinh(2(^). 



sinh(2(^). 



(2.1) 
(2.2) 



All quantities in the physical frame are denoted with a tilde, and any quantity without a tilde is in the Einstein frame. 
The total action of TeVeS, 5, contains contributions from the three dynamical fields and a matter contribution (see 
Q for details). These include two positive dimensionless parameters, k and K , which are the coupling parameters for 
the scalar and vector fields respectively. There also exists a dimensionless free function _F, a constant length scale 
and a spacetime dependent Lagrange multiplier, A. 

By varying the total action, S, with respect to g^'^ , one can obtain the field equations for the tensor field 



Guv — SttG 



(2.3) 



where T^i, is the energy-momentum tensor in the physical frame, T^(^^U^-^ = Ta^JAv + Ta^hi^ and G^j/ is the Einstein 
tensor in the Einstein frame. Conservation of energy-momentum is therefore given in the physical frame as V/^T'"' — 0. 
The other sources in Ea. (|2.3p are given by 



Ga2 



F{kGa'^)gf,^ -W^tp^a (l^{f,V,u) - ]^U^Vj39iJ^i^ 



(2.4) 
(2.5) 



where ^q,^] = Wa,/3 —l^p,a- Similarly, by varying S with respect to lA^ and Lp, one obtains the field equations for the 
vector and scalar fields; 



+ + ^TTGa^U'^p^f3g°'''ip^^ = SvrG (f - e"^'^) g^'^U'^f^p, 
[^(M2/i^>,^^,,)/i«/'^,„] = kG [g'^^ + {1 + U^U"] f^p, 



(2.6) 
(2.7) 



where h"^ 
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-W^U^ and ^{x) is a function defined by 2/iF(/i) -f ^^dF{fj,)/dii — —2x. With this function /i, the 



nondynamical scalar field a is determined by 



kGa^ = fi{ke^h'"^(p^o.p.,0)- 



(2., 
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Therefore, the field equations of TeVeS are Eqs . (|2.3p and (|2.6p - (|2.8p . It has been shown in the strong-field limit 
that /J, = 1 is an excellent approximation [ill, [12] ■ On cosmological scales this is not a good choice .Gj] , however in this 
article we only consider regions not too far from neutron stars, and we therefore set /i = 1. This implies from Eq. 
(|2.8p that cr^ — \/{kG). Moreover, whilst the functional form of F is not predicted by the theory, one can show that 
when /i = 1, the contribution of F to the field equations vanishes [l,[ill[i2|- Therefore, our results are independent 
of this function and we drop it from the remaining discussion. 



B. TOV in TeVeS 



First, the Tolman-Oppcnheimer-VolkofF (TOV) equations in TeVeS are derived by Lasky, Sotani & Giannios [l^ 
Here we make an brief description of TOV equations. A static, spherically symmetric metric can be expressed as 



9af3 



dx^dx' 



(2.9) 



where d^"^ = dO"^ + s\r? 6d(j)'^ and e^^ = 1 — 2ni{r) jr. In general, the vector field for a static, spherically symmetric 
spacetime can be described as = (W*,W",0,0), where and W are functions of r. Giannios [ll| showed that in 
vacuum, the parameterized post-Newtonian (PPN) coefficients for a spherically symmetric, static spacetime with a 
non-zero can violate observational restrictions. In this article, we therefore only consider the case where W = 0. 
In this case, the vector field can be fully determined from the normalization condition, and is found to be = 
{e'"/'^, 0,0,0). Moreover, one can show that the vector field equation (|2.6[) is now trivially satisfied. With this vector 
field, the physical metric is 



ds'^ = gap 



dx'^dx^ 



(2.10) 



and the fluid four- velocity is = e'^Z//^. We further assume the stellar matter content to be a perfect fluid, i.e., 
Tp,i, = {p + P)u^Ui, + PQiivi from which one can show that the full system of equations with k ^ and K ^ reduces 
to 



1-^ 
2 
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-2ip 



Kr 
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27rr^ 
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{2ij + j,'), 



(2.11) 
(2.12) 

(2.13) 
(2.14) 

(2.15) 



where a prime denotes a derivative with respect to r. (See [l3 | for the derivation of these equations and for a discussion 
with k = and/or with K = 0.) This system of equations is closed with the addition of an equation of state (EOS). 
The stellar radius in physical frame, R, is determined by i? = e~'^^'^'^rs, where is the position of the stellar surface 
defined as the point where P = 0. Note that on exterior region the scalar field still exists although there is no fluid. 

We integrate the above equations from the center, r = 0, to the stellar surface, r = Vs- Moreover, the interior 
boundary conditions are given by P{0) = Pq, p{0) = po, v{0) = vq, </?(0) = (po, • 0(0) = 0, and m(0) = 0, which are 
determined by Taylor series expansions of the above equations near r = (see 1J| for details). The exact values 
for vq and ipo are determined by matching the functions ^{r) and f{r) to their asymptotic behavior, which is found 
by performing a coordinate transformation on (|2.10p to bring it into an asymptotically flat form. We define new 



coordinates t = te'^" and 



re 



= , where (p^ denotes the cosmological value of the scalar field. Then performing 



an asymptotic expansion of all the equations and dropping the hats on the new coordinates for simplicity in the 
expressions implies 



2Madm , Z' ^ 
gtt = -IH h O 



9rr 1 



2Madm 



O 



(2.16) 
(2.17) 
(2.18) 
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Here, -Madm is the total Arnowitt-Deser-Misner (ADM) mass given by 

Mabm = (m^ + (2.19) 

where TOoq is the mass function evaluated at radial infinity. Also, M^p is the scalar mass @, which is constant outside 
the star and is defined everywhere as 

^Att J + 3P) e'^-'+^^^^-^'^dr. (2.20) 

We adopt the same EOS as in jSjl, which are polytropic ones derived by fitting functions to tabulated data of realistic 
EOS known as EOS A and EOS II. The maximum masses of neutron stars in GR are M — 1.65Af0 with R = 8.9 km 
for EOS A and M = 1.95Mq with R = 10.9 km for EOS II. That is, EOS A is considered soft, while EOS II is an 
intermediate EOS. 

When it comes to the study of the structure of neutron stars, TeVeS introduces three new parameters, k, K, and 
(/5c, with respect to GR. Since the value of k is tightly constrained by both cosmological models and also planetary 
motions in the outer solar system Q, we accordingly set k = 0.03 for the remainder of the article. With respect 
to the value of <pc, Lasky, Sotani & Giannios showed that can have a minimum value of around 0.001, based on 
causality issues inside the neutron star [l3|. Therefore, for this article we use tp^ = 0.003. Details of neutron star 
models where these parameters are allowed to vary are given in Lasky, Sotani & Giannios (l3 | , in which they showed 
that the dependences on k and on are minimal for neutron star models. While, restrictions on K are less severe, 
and have not been discussed in great detail in the literature. In this article we consider the range < X < 2, because 
for i^T > 2 one can show that the pressure diverges from the stellar center outward, and therefore stellar models are 
not possible [l3| (Sagi & Bekenstein [l3| also showed that physical black hole solutions are only valid for K <2). 



C. Neutron Star Models in TeVeS 



Fig. [T] shows the ADM mass as a function of the central density of neutron stars. Different lines correspond to 
different values of K, whose values are indicated. Additionally, we plot the stellar model for GR with the solid line. 
Note that for spherically symmetric neutron stars, stellar models for the region dMADu/dpo < could be unstable. 
From this figure we can see that, although the central density giving the maximum mass is almost independent from 
the value of K, the corresponding maximum mass depends strongly on the parameter i.e., for larger values of K 
the maximum mass becomes smaller. For example, with EOS A for K = 0.5 the maximum mass is 18 % smaller than 
that of GR. 




FIG. 1: Relation between the mass and central density of neutron stars in GR and in TeVeS with k = 0.03, ifc = 0.003, where 
the left and right panels are corresponding to the stellar properties given by EOS A and EOS II, respectively. In the figure 
the solid line denotes the case of GR while the other lines are corresponding to the stellar models with different values of K in 
TeVeS. 

Fig. [2] shows the relation between the ADM mass and stellar radius with different values of K and also for the GR 
case. In general, one requires a softer EOS near the stellar surface, which implies the stellar radius becomes larger, but 
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with high central density for typical neutron star, the stellar models are almost independent from the consideration 
of softer EOS near the stellar surface. This figure further implies that there exists a minimum radius, which usually 
corresponds to the maximum mass [l^. Considering this minimum radius, it can be seen in Fig. [2] that neutron stars 
in TeVeS are smaller than in the GR case. For example, we can see that the minimum radius for a star with EOS A 
in TeVeS with K = 0.5 is 7.7 km, whereas for GR it is 8.9 km, which is a 13.5% difference. 
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FIG. 2: Mass-radius relation for neutron stars in GR and in TeVeS with k = 0.03, v^c = 0.003, where the left and right panels 
are corresponding to the stellar properties given by EOS A and EOS II, respectively. 



III. PERTURBATION EQUATIONS IN THE COWLING APPROXIMATION 

In this section we derive the perturbation equations for nonradial oscillations of spherically symmetric neutron stars 
in TeVeS. For simplicity, we adopt the Cowling approximation, in which the fluid is perturbed on a fixed background. 
That is, the perturbations of the spacetime, vector field and scalar field are frozen, i.e., (5g^i, = 0, Sli^ — and 
= 0. It notes that with the Cowling approximation we can study only the oscillation modes related to the fluid 
perturbations, such as /, p and g modes, while it is impossible to study the other emissions of scalar waves, vector 
waves and gravitational waves connected to the oscillation of spacetime. Additionally, we should notice that the 
Cowling approximation in GR is typically very good for axial type of oscillations while for polar type of oscillations 
the error for typical rclativistic stellar models could become less than 20 per cent for / modes and around 10 per cent 
for -p modes [2J] . 

With Cowling approximation, the perturbed energy- momentum tensor, (JT^", is given as 

5f '"^ = {8~p + u^u"" + + p) + u^iu'') + bP-g^". (3.1) 

Introducing the Lagrangian displacement vector, the perturbed variables in bT^^ such as ^u^, 8p and 5P, can be 
described explicitly. The Lagrangian displacement vector for the fluid perturbations are 

r = (r, l\ 1^) = (W^, -Vde. -Fsin"2 m^) -^F,™, (3.2) 

where W and V are functions of t and r. Then the perturbations of 4-velocity, can be written as 

bu^ = (o,VF,-l/a0,-ysin-2 0a^j e-'^-'^/^^r,™, (3.3) 

where dots on the variables denote the partial derivative with respect to t. On the other hand, using the first law 
of thermodynamics, we can get the following relation between the adiabatic changes of the density and the baryon 
number density; 



(3.4) 
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where n denotes the baryon number density. So if we use the relationship between the Lagrangian perturbation Ap 
and Eulerian perturbation dp such as 



A/5 ~ (5p + C^rP, 



we can express the Eulerian density variation as 



~,An p'W^^ 
5p^{p + P)— - ^Y£„ 



Additionally, with the definition of the adiabatic constant 

/ ainp' 
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(9 Inn 



hAP 
PAn 



we can derive the Eulerian variation of the pressure; 



~ All P'W 
6P = jP- — 

n 

It notes that with Eqs. (|3.4p and (|3.7p we can get the useful expression for 7 as 

p + P f dP\ 



P 



dp 



(3.5) 
(3.6) 

(3.7) 

(3.8) 
(3.9) 



Finally, the Lagrangian variation of the baryon number density, which comes on the expressions of Sp and 5P, is 
determined by the relation as 



An _ _^{3)Xk _ ^9 



(3.10) 



~ (3) 

where and g denote the covariant derivative in a 3-dimension with metric 5^1/ and the determinant of gf^^, 

respectively. In this article, since we assume the Cowling approximation the second term is neglected. Then the 
Lagrangian variation of the baryon number density can be written as 



An 

n 



^' + ^ (C - Gf') W + l{£+ 1)V 



1 



tY, 



im ■ 



(3.11) 



Finally we can get the equations describing the fluid perturbations by taking a variation of the energy-momentum 
conservation law, W^T^" — 0. With Cowling approximation, this equation becomes V^ST'^'^ = 0. The explicit forms 
with p = r,6 are 



dr 



jP 



1 1 P'W 



p' + P' 



p + P 



W + - (C' - 6(^') W + l{t+ \)V \ + 



P'W 



w 



0, 



IP 



w + - iC ~ 6^') w + ei£ + i)v } + 



P'W 



= 0, 



(3.12) 



(3.13) 



where we use the relation that Sp/SP — p' / P' and Eq. (j2.13p . By assuming a harmonic dependence on time, the 
perturbative variables will be written as W{t,r) — W^(r)e*'^* and V{t,r) — F(r)e*'^*. To make the above equations 
simpler, by calculating the combination of the form (j (|3.13[) /(jr — (I3.12p and substituting Eq. (I3.13|) again, we can get 



V = (V + iy')V - eS 



W 



(3.14) 



Thus, from Eqs. p.l3p and (|3.14p . we can obtain the following simple equation system for the perturbations of fluid; 



W' = ^ 
dP 



W 



+ -{6ip' ~C)w-e{£ + i)v, 



V = [A^' + v')V -e^ — 



(3.15) 
(3.16) 
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With the appropriate boundary conditions at the center and stellar surface, the above equation system constitutes 
an eigenvalue problem for the parameter to. One can find the behavior of W and V near the stellar center as 
W{r) = Br^^^ + 0{r^^^) and V{r) = —Br^/£ + 0(r^+^), where B is an arbitrary constant, while the boundary 
condition at the stellar surface is the vanishing the Lagrangian perturbation of the pressure, i.e., AP = 0. Since the 
Lagrangian perturbation of the pressure is described by AP = 'yPAh/h, with a help of Eq. p.lSp we can get the 
boundary condition at the stellar surface as 

2oj^r^e^^f-''V + {2(p' + i'')W ^0. (3.17) 



IV. OSCILLATION SPECTRA 



With respect to the neutron star models shown in Sec. m in this section we examine the stellar oscillations. 
Especially, we focus on the stellar models whose central density is in the range from po = lO^"^ g/cm^ up to the value 
given the maximum ADM mass. The stellar parameters with maximum ADM mass are summarized in Tables HI and ITll 
where tpo and z are central value of ip and surface redshift, respectively. In general, the oscillation spectrum is directly 
related to the stellar parameter, such as mass, radius and EOS, but the frequencies of fundamental oscillation modes, 
i.e., / modes, can be connected to the stellar average density, {Mabm/R^Y^^- This reason is physically explained by 
considering the relation between the sound speed and the time that fluid perturbation needs to propagate across the 
star. Actually, for the stellar models in GR, Andersson & Kokkotas found the empirical formula for the frequency of 
/ mode as a function of stellar average density [1^ . The / mode frequencies for the stellar models in GR constructed 
with almost all EOS are subject to this empirical formula. While, the frequencies of / mode for the stellar models in 
TeVeS with above two EOS are shown in Fig. [31 The deviation from GR is clearly cognized for typical neutron stars 
and depending on the value of parameter the frequencies become around 20 % larger than those expected in a 
general relativistic neutron star. This can be an observable effect and one might distinguish the gravitational theory 
in strong gravitational field by using the observations of gravitational waves. 



TABLE I: Stellar parameters for models with EOS A and with maximum ADM mass, where we choose that k — 0.03 and 
ipc = 0.003. 



K 


Madm/Mq 


po [g/cm^] 


R [km] 


fo 




Madm/R 


z 


0.2 


1.56 


3.54 X lO^'^ 


8.36 


8.10 X 10" 


-4 


0.271 


0.479 


0.5 


1.42 


3.30 X lO^'^ 


7.60 


1.26 X 10" 


-3 


0.265 


0.459 


1.0 


1.15 


3.00 X 10^^ 


6.16 


1.91 X 10" 


-3 


0.255 


0.433 


1.5 


0.81 


2.81 X 10^^ 


4.32 


2.48 X 10" 


-3 


0.248 


0.414 


1.9 


0.36 


2.70 X 10" 


1.92 


2.90 X 10" 


-3 


0.242 


0.401 



TABLE II: Stellar parameters for models with EOS II and with maximum ADM mass, where we choose that k = 0.03 and 

Lpa = 0.003. 



K 


Madm/Mq 


Po [g/cm^] 


R [km] 






Madm/R 


z 


0.2 


1.84 


2.40 X 10" 


10.32 


9.38 X 10" 


-4 


0.259 


0.441 


0.5 


1.67 


2.25 X 10" 


9.39 


1.35 X 10" 


-3 


0.253 


0.425 


1.0 


1.36 


2.06 X 10" 


7.63 


1.96 X 10" 


-3 


0.245 


0.403 


1.5 


0.96 


1.92 X 10" 


5.36 


2.51 X 10" 


-3 


0.238 


0.385 


1.9 


0.43 


1.80 X 10" 


2.40 


2.91 X 10" 


-3 


0.231 


0.370 



The possibility to probe the gravitational theory by using observations of gravitational waves can be also seen in 
Fig. m where we plot the normalized frequencies of / and pi modes as functions of ADM mass. In the figures, the 
solid line denotes the frequency in GR while the other broken lines are corresponding to those in TeVeS with several 
values of K. One can easily observe that the frequencies expected in TeVeS are quite different from those in GR. Since 
this distinction results from the difference of gravitational theory, which creates due to the presence of scalar field, 
observing more than one mode of gravitational wave could tell us the existence of the scalar field. This statement 
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FIG. 3: The frequency of / mode as a function of the stellar average density, {Madm/ R^)^^^ , where // is defined as // = 
u}f /{2ti). The solid line corresponds to the frequency in GR, while the other broken lines are corresponding to the frequencies 
in TeVeS with several values of K. Notice that the unit of average density is [1/km] in the geometrical unit, where c = G = 1. 



might become more obvious by seeing the dependence of frequencies of gravitational waves on the parameter K . In 
Fig. O we plot the normalized frequencies of the first four modes, i.e., /, pi, p2 and ps, as functions of parameter 
K , where the ADM masses are fixed to be 1.4M0. The allowed maximum values of K to produce the stellar models 
with Madm = 1.4M0 are K = 0.54 for EOS A and K = 0.94 for EOS II. From these figures, we can see that the 
qualitative dependences of frequency on the value of K are independent from EOS and kinds of eigen-mode. That 
is, the normalized frequencies of fluid modes are decreasing as the value of K becomes large. On the other hand, it 
is also found that the higher overtone is quantitatively more sensitive against the value of K than the lower modes. 
This point can be seen in Table lllli where we summarize the ratio of difference between the frequencies for the stellar 
models with K — 0.05 and with the allowed maximum values of K. In other words, the frequencies of higher overtone 
is more helpful to distinguish TeVeS from GR via the gravitational wave observations. Anyway, through Figs. [4] and 
[5l we can find that with a help of observation of stellar mass, it is possible to probe the gravitational theory in the 
strong-field regime by using observations of gravitational waves. 



TABLE III: The relative frequency change of each eigen-mode in Fig. [5] defined as ((x;/f=o.05 — L^K=max) /l^k =0.05, where 
ujk^o.o5 and uoK^max denote the frequencies for the stellar models with K = 0.05 and with the allowed maximum values of K, 
respectively. 



mode 


EOS A 


EOS II 


/ 


10.28% 


12.89% 


pi 


15.21% 


18.18% 


P2 


17.88% 


20.73% 


PS 


19.23% 


21.89% 



V. CONCLUSION 



In this article, to examine the effect of the Tensor- Vector- Scalar (TeVeS) Theory on the oscillation spectra of neutron 
stars, we have derived the perturbation equations of neutron stars in TeVeS and calculated their eigen-frequencies. 
Depending on the parameter of TeVeS, the frequencies of fundamental oscillation could be off the well-known empirical 
formula in GR and they become lager than those expected in GR. We can also see the deviation from GR in the 
frequencies of higher overtones and they have stronger dependence on the parameter K than the lower oscillation 
modes. Since these imprints of TeVeS come from the presence of scalar field, by using the observations of gravitational 
waves associated with the stellar oscillations, it will be possible not only to distinguish the gravitational theory in the 
strong-field regime, but also to probe the existence of the scalar field. 

For simplicity, we assumed the Cowling approximation in this study, which restricts our examination to only 
stellar oscillations. This means that we should do more detailed study including the metric, vector and scalar fields 
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FIG. 4: The normalized frequencies of first two fiuid modes are plotted as functions of the ADM mass, where upper and lower 
panels are corresponding to / and pi modes, respectively. In these figures, the solid line corresponds to the frequency in GR, 
while the other broken lines are corresponding to the frequencies in TeVeS with several values of K. 



perturbations. Via these oscillations, we could obtain the additional information in the gravitational spectrum, and 
combining those with results shown in this article would provide more accurate constrains on the gravitational theory 
in the strong-field regime. Furthermore the introduction of the stellar magnetic effect might be also important. For 
example, recent observation of quasi-periodic oscillation in the giant flares are believed to be related to the oscillations 
of strong magnetized neutron stars [2^, [13, HI] . Considering the magnetic effects, one might be able to get the further 
constraint in the theory. 



Acknowledgments 

We thank K.D. Kokkotas for valuable comments. This work was supported via the Transregio 7 "Gravitational 
Wave Astronomy" financed by the Deutsche Forschungsgemeinschaft DFG (German Research Foundation). 



[1] D. Psaltis, preprint (0806.1531 [astro-ph]). 

[2] S. DeDeo and D. Psahis, Phys. Rev. Lett. 90, 141101 (2003). 

[3] H. Sotani and K.D. Kokkotas, Phys. Rev. D 70, 084026 (2004); 71, 124038 (2005). 
[4] T. Damour and G. Esposito-Farese, Classical Quantum Gravity 9, 2093 (1992). 
[5] CM. Will, Living Rev. Relativity 4, (2001). 
[6] J.D. Bekenstein, Phys. Rev. D 70, 083509 (2004). 
[7] M. Milgrom, Astrophys. J. 270, 365 (1983). 
[8] D.M. Chen and H.S. Zhao, Astrophys. J. 650, L9 (2006). 
[9] C. Skordis et al., Phys. Rev. Lett. 96, 011301 (2006). 
[10] S. Dodelson and M. Liguori, Phys. Rev. Lett. 97, 231301 (2006). 



10 




FIG. 5: For the stellar models with Madm = IAMq, the normalized eigenvalues ui of the first few modes (/, pi, P2 and pa) 
are shown as functions of parameter K with EOS A (solid lines) and EOS II (broken line). The values of k and ipc are fixed as 
k — 0.03 and ip^ = 0.003, respectively. 



[11] D. Giannios, Phys. Rev. D 71, 103511 (2005). 

[12] E. Sagi and D. Bekenstein, Phys. Rev. D 77, 024010 (2008). 

[13] C.R. Contaldi, T. Wiseman, and B. Withers, Phys. Rev. D 78, 044034 (2008). 

[14] P.D. Lasky, H. Sotani, and D. Giannios, Phys. Rev. D 78, 104019 (2008). 

[15] N. Andersson and K.D. Kokkotas, Phys. Rev. Lett. 77, 4134 (1996). 

[16] H. Sotani, K. Tominaga, and K.I. Maeda, Phys. Rev. D 65, 024010 (2001). 

[17] H. Sotani and T. Harada, Phys. Rev. D 68, 024019 (2003); H. Sotani, K. Kohri, and T. Harada, ibtd 69, 084008 (2004). 

[18] A. Stavridis, A. Passamonti, and K.D. Kokkotas, Phys. Rev. D 75, 064019 (2007). 

[19] M. Vavoulidis, K.D. Kokkotas, and A. Stavridis, Mon. Not. R. Astron Soc. 384, 1711 (2008). 

[20] H. Sotani and M. Saijo, Phys. Rev. D 74, 024001 (2006). 

[21] H. Sotani, S. Yoshida, and K.D. Kokkotas, Phys. Rev. D 75, 084015 (2007). 

[22] C.M. Will, TheoryandExperimentinGravitationalPhysics (Cambridge University Press, Cambridge, 1993). 

[23] J.M. Lattimer and M. Prakash, Astrophys. J. 550, 426 (2001). 

[24] S. Yoshida and Y. Kojima, Mon. Not. R. Astron Soc. 289, 117 (1997). 

[25] N. Andersson and K.D. Kokkotas, Mon. Not. R. Astron Soc. 299, 1059 (1998). 

[26] H. Sotani, K.D. Kokkotas, and N. Stergioulas, Mon. Not. R. Astron Soc. 375, 261 (2007); 385, L5 (2008). 

[27] H. Sotani, A. Colaiuda, and K.D. Kokkotas, Mon. Not. R. Astron Soc. 385, 2161 (2008). 

[28] H. Sotani and K.D. Kokkotas, accepted in Mon. Not. R. Astron Soc, preprint (0902.1490 [astro-ph]). 



